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(N Abstract 

Future galactic supernovae will provide an extremely long baseline for studying the properties and 
interactions of neutrinos. In this paper, we discuss the possibility of using such an event to constrain 
(or discover) the effects of exotic physics in scenarios that are not currently constrained and are not 
{N) accessible with reactor or solar neutrino experiments. In particular, we focus on the cases of neutrino 

decay and quantum decoherence. We calculate the expected signal from a core-collapse supernova 
in both current and future water Cerenkov, scintillating, and liquid argon detectors, and find that 
such observations will be capable of distinguishing between many of these scenarios. Additionally, 
future detectors will be capable of making strong, model-independent conclusions by examining events 
Ph associated with a galactic supernova's neutronization burst. 
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I. INTRODUCTION 



In many scenarios of physics beyond the standard model, the associated signatures are 
observable only in high energy interactions, motivating the use of particle accelerators such as 
the Large Hadron Collider and the Tevatron. In other cases, however, the effects of new physics 
can appear over very long periods of time, or in very long baseline experiments. The discovery 
of neutrino flavor oscillations provide such an example. Neutrino oscillations have only been 
observed among neutrinos propagating between the Sun and Earth and over shorter distances, 
however, leaving open the possibility that other new phenomena might become observable 
through the study of neutrinos traveling over much greater distances [H [2]. 

A wide variety of distant astrophysical objects are predicted to potentially produce observ- 
able fluxes of neutrinos, including active galactic nuclei, gamma ray bursts, microquasars, X-ray 
binaries, and supernovae [3]. Among these prospective sources, however, only neutrinos from 
supernovae have been observed. In particular, a total of 24 (likely anti-electron) neutrinos were 
detected by three different detectors in the seconds following supernova 1987A. At a distance of 
approximately 50 kiloparsecs from Earth, SN 1987A provided an opportunity to observe neu- 
trinos over a baseline that is roughly 10 10 times longer than that traveled by solar neutrinos. 
When the next nearby supernova takes place, modern detectors will observe far more events 
and collect a much greater range of information than was possible in 1987. 

In this article, we discuss the possibility of using neutrinos from a galactic supernova as 
a probe of exotic physics. In particular, we focus on two examples of phenomena that are 
observable only over long baselines or timescale: neutrino decay and quantum decoherence. The 
remainder of this article is structured as follows. In Sec. [TTJ, we summarize the exotic physics 
scenarios considered in our study. In Sec. Ill, we discuss the predicted neutrino spectrum from 
supernova and in Sec. IV we describe the experiments used to detect them. In Sec. [Vj we 
describe how neutrino decay or quantum decoherence would impact the neutrinos observed 
from a galactic supernova and discuss how such observations could be used to constrain or 
identify such phenomena. In Sec. VI we discuss the role that neutrinos from a supernova's 



neutronization burst can play in identifying or constraining new physics. Finally, in Sec. VII 
we summarize our results and conclusions. 



II. EXOTIC PHYSIC SCENARIOS 

A. Neutrino Decay 

In this study, we focus on two body decays of the form 

Vi — > Vj + X, (1) 

where z/j and Uj denote either neutrino or anti-neutrino mass eigenstates and X represents 
a very light or massless particle (such as a majoran) that is not detected (for related model 
building, see Ref. [1]). In contrast to radiative decays [5], decays of this form are viable over 
supernova-scale baselines. The strongest constraints on decays of this type are on the order of 
r/m > 10~ 4 s/eV, and result from solar neutrino data [6] and from searches for the appearance 
of v e at KamLAND [7j. In some cases, constraints can also be placed from the observations of 
SN 1987A. 
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The 24 neutrino events associated with SN 1987A marked the first (and thus far only) 
detection of neutrinos from a supernova. Although the modest number of detected events limit 
the utility of any extensive statistical analysis, some general features have been confirmed. 
Specifically, the number of events detected from SN 1987A agree well with expectations. From 
this, we can conclude that z/ e 's do not significantly decay, at least over timescales of less than 
~ 1.7 x 10 5 years, which limits the corresponding rest-frame lifetime to r D< , > 0.7 days x 

/ m Vf , \ ( 10MeV\ 
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The constraints from SN 1987A are most relevant for the case of an inverted mass hierarchy, 
in which u 3 is the lightest mass eigenstate and, due to the small value of #13, has a very small 
projection onto u e . Even if # 13 is taken to have the largest experimentally allowed value, the 
decay scenario where both v\ — > v 3 + X and v% — > V3 + X leads to about a factor of 10 
suppression in the number of expected events. As this suppression was not observed in SN 
1987A, we can conclude that the timescale for this process (for E u ~ 10 MeV) must also be at 
least ~ 1.7 x 10 5 years. The baseline to the Large Magellanic Cloud (where 1987A occurred) is 
extremely long at approximately 50 kpc. As this study is focused on galactic supernova with a 
baseline of roughly 10 kpc, this scenario of neutrino decay can be excluded from our analysis. 
Within the context of the normal hierarchy, however, SN 1987A does not strongly constrain 
the possibility of decays of one neutrino mass eigenstate to another. 

In the decay of a neutrino, helicity can either be conserved or flipped. If the masses are 
quasi-degenerate (m^ — rrij <C mi,mj), then the conserved helicity decay is strongly favored. 
In this case, the daughter neutrino receives nearly all of the parent's energy, and the overall 
spectrum is left nearly unchanged. If the masses are instead hierarchical (rrii — rrij ~ rrii), 
then a larger fraction of the parent energy can be lost to the X particle, leading to a degraded 
spectrum of the daughter species relative to that of the parent. 

The details of the resulting spectrum further depend on whether the neutrinos are Majorana 
or Dirac particles. If the neutrinos are Majorana particles, then the daughters with flipped 
helicity will be detectable as a corresponding anti-neutrino. With the parent (daughter) energy 
labeled as E p (E d ), the daughter spectrum will be scaled down by a factor of (E p — E d )/E^ 
relative to that of the parent [8] . The daughters with conserved helicity will also have a spectrum 
found by scaling that of the parent, but this time by a factor of E d /E p [8]. If the neutrinos 
are Dirac particles, then the daughter neutrinos with flipped helicity will be sterile and thus 
undetectable. The daughters with conserved helicity will again have a spectrum produced by 
scaling the parent spectrum by a factor of E d j 1 E p . Thus in such a decay scenario, not only will 
the energy spectrum be degraded, but also the total number of neutrinos/anti-neutrinos will 
be suppressed. 

B. Quantum Decoherence 

Within the context of standard quantum mechanics, a pure state will indefinitely remain in 
that state, never oscillating into a superposition or mixture of states. It has been suggested, 
however, that quantum fluctuations of the gravitational field could potentially alter this conclu- 
sion. Short lived microscopic black holes, for example, could appear as part of the space time 
foam, leading to the loss of quantum information as particles move across their event horizons, 
and transforming pure states into a mixed ones [9] . 

Among neutrinos, one might hope to observe quantum decoherence through the loss of flavor 
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information taking place over very long baselines [2J [10]. The prediction for this phenomenon 
is the evolution of neutrino flavors toward a ratio of v e : : v T — 1/3 : 1/3 : 1/3, for both 
neutrinos and anti-neutrinos, regardless of their initial conditions. If this evolution takes place 
over kiloparsec distance scales, it may be possible to observe with the next galactic supernova. 

III. NEUTRINOS FROM SUPERNOVAE 

Over a period of approximately 10 seconds, a supernova produces approximately as many 
neutrinos as the Sun will generate over its entire life. Approximately 99% of the gravitational 
energy associated with a stellar core-collapse supernova is released in the form neutrinos. 

Despite considerable effort spanning several decades in the modeling of core-collapse super- 
novae, there are still many aspects of the explosion that are not well understood. In 1998, the 
Livermoore group produced the most recent simulation that obtained an explosion and followed 
the evolution through 18 s [11]. The model predicts strongly hierarchical time-integrated, av- 
erage energies: (E Ue )?al2 MeV, (E De )^15 MeV, and (E Ux )^2A MeV, where u x denotes u^, 
u T , u^, and v T . This strong hierarchy is useful for distinguishing between different oscillation 
and decay scenarios, as there are several important neutrino interactions involving v e and v e 
with relatively high energy thresholds that would not produce many events without significant 
mixing from the harder v x spectra. The Livermoore group also found luminosity and total 
energy equipartition between the flavors throughout the evolution. 

The time integrated neutrino spectrum is often written as a "pinched" Fermi-Dirac distri- 
bution: 

dNt, = NEl 

dE v T 4 r e ^-v + 1 y u 

with E v the energy of the neutrinos, N an overall normalization constant, T the neutrino 
temperature, and rj is a chemical potential for the neutrinos. For this spectrum, rj is a measure 
of how pinched the low energy and high energy tails would be when compared to a Fermi-Dirac 
spectrum with the same average energy. The parameters that we use in our model are chosen 
to match well with the spectra of the Livermore model, with T Ue =3.5 MeV, T Pe = 4 MeV, 
T Ux = 7 MeV, i] Ue = 2.08, r] Pe = 2.5, and r\ Vx = [IT] . These spectra are shown in Fig. [I] 

Although the Livermoore group did obtain an explosion with their simulation, the model 
neglected neutrino processes that are now thought to be important. More recent simulations 
which account for these processes fail to produce an explosion and only follow the supernova 
evolution for a maximum of ~1 s (see Table III of Ref. [J2] for references to many of these 
models). As discussed in Ref. [T2], luminosity equipartition between v e and v e seems to hold 
for most models, but u x luminosity can vary by a factor of 2 in either direction depending on 
the model and evolutionary phase of the supernova. Additionally, most of these models produce 
a somewhat less pronounced energy hierarchy than that found in the Livermoore model. For 
example, the Garching Group reports (E Ue ) « 12 MeV, (E De ) w 15 MeV, and (E Ux ) w 18 MeV 

Em 
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FIG. 1: The spectrum of neutrinos from a supernova, prior to oscillations, as described by Eq. [2] 

IV. THE DETECTION OF SUPERNOVA NEUTRINOS 

In this section, we summarize the primary detection techniques used to potentially detect 
the neutrinos from a nearby supernova. 

A. Water Cerenkov Detectors 

When charged particles such as electrons travel faster than the local speed of light, they 
emit radiation which can be detected. Currently, the largest water Cerenkov detector is Super 
Kamiokande. The signal from an individual supernova is a short lived event (~10 s), which 
permits higher background rates, allowing Super Kamiokande to use its full fiducial mass of 
32 kilo-tons (kton), significantly larger than any other existing neutrino detector, although 
detectors at the 500 kton scale have been proposed (UNO [15] . MEMPHYS [15], and Hyper- 
Kamiokande [16]). The interactions relevant for the observation of supernova neutrinos by 
water Cerenkov detectors include: 

1. Inverse Beta Decay: v e + p — > n + e + 

For any target material with a large number of free protons (those not bound in a nucleus), 
inverse beta decay events will typically dominate the signal. This is due in large part to the 
relatively large cross section and the low energy threshold (~1.8 MeV). This reaction has been 
studied extensively both experimentally and theoretically. We use the cross section from Ref. 
[17] which includes radiative corrections and an exact treatment of the kinematics. 

A proposal has been made to dissolve approximately 0.2% gadolinium within the water of 
Super Kamiokande with the purpose of tagging the neutrons produced by inverse beta decay 
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[18J. These neutrons quickly thermalize and are captured on the gadolinium nucleus with a 
90% efficiency, producing a cascade of detectable photons [T9] . 

2. Electron Elastic Scattering: v + e~ — > v + e~ 

Although all flavors of neutrinos can undergo the process of elastic scattering with electrons, 
the cross sections vary for the different flavors. This can be traced to the fact that for electron 
neutrinos both charged current and neutral current processes contribute to the cross sections. 
Supernova neutrinos have energies far below the mass of the muon, so only the neutral current 
process contributes to the cross section for v x . The approximate ratios of the total cross sections 
are a Ve : a De : a Ux :a- Ux -l: 0.42 : 0.16 : 0.14 [2U]. 

Although there are copious numbers of electrons in any target material, the cross section is 
typically much smaller than for interactions with protons or neutrons, causing the total event 
rate to be dominated by nuclear scattering. Electron scattering events are strongly forward 
peaked (due to the small electron mass), however, providing a means by which to separate these 
events from the large rate of inverse beta decay. If the two signals can indeed be separated 
then this process could provide an opportunity to measure the u e component of the neutrino 
spectrum. 

3. Scattering With Oxygen Nuclei: u e + w O -> 16 N+e + ; u e + 16 -> 16 F+e~, ^+0 -> v+O* 

Supernova neutrinos can also interact with oxygen nuclei in water Cerenkov detectors. 
Both v e and u e undergo charge current interactions with oxygen nuclei and have thresholds 
of 15.4 MeV and 11.4 MeV, respectively. The signals from the resulting positrons and elec- 
trons are indistinguishable from each other and from inverse beta decay positrons in Super 
Kamiokande, so these reactions can not be uniquely identified. Additionally, the high thresh- 
olds mean that unless the harder v x spectra mixes significantly with the v e and v e (through 
oscillations or other exotic physics), the oxygen charged current signal will be quite low. 

All neutrino flavors can undergo neutral current reactions with oxygen nuclei, which can 
yield detectable /3 and 7 particles in the subsequent decays of the excited nucleus. The most 
interesting of these decay chains involve excited states of 15 and 15 N . These states are only 
able to decay through the emission of a 7 with an energy in the range of 5 MeV to 10 MeV and 
are detectable as peaks above the smooth background (see Figure 2 of Ref. [2TJ). This oxygen 
reaction was suggested as a supernova neutrino signal in Ref. [22], where they estimated 300 
oxygen events over a smooth background of roughly 270 inverse beta decay events in this energy 
window. 

B. Scintillating Detectors 

Scintillating detectors generally contain an organic oil that is composed primarily of carbon 
and hydrogen. Currently, KamLAND and LVD are the largest scintillating detectors with a 
mass on the order of 1 kton [231 124"] . Although these detectors are much smaller than Super 
Kamiokande, scintillators typically have much lower energy thresholds and higher light output 
allowing them to compete with Super Kamiokande in certain arenas. Proposals for larger 
scintillators include LENA (50 kton) |25j and Hanohano (10 kton) [26] . 
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As scintillators contain many free protons, inverse beta decay provides a strong supernova 
signal in these detectors. These events can be tagged in scintillating detectors as the positron 
produced annihilates with an electron producing a pair of photons that are above the energy 
detection threshold. Additionally, there will be a delayed coincidence signal from a 2.2 MeV 
photon produced when the neutron captures onto a proton producing deuterium. As with 
water Cerenkov detectors, electron elastic scattering provides a much smaller signal. We will 
now briefly describe the other important interactions for scintillators. 

1. Scattering With Carbon Nuclei: u e + 12 C ->■ 12 B+e+, z^ e + 12 C -)• 12 N+e", v+C ->■ v+C* 

At supernova neutrino energies, both v e and v e can undergo charged current reactions. 
The v e reaction has a threshold of 17.3 MeV and produces 12 N which subsequently decays by 
emitting a e + with a mean lifetime of 11.0 ms [23J. The emitted positron annihilates producing 
a gamma ray pair that give a delayed coincidence signal for tagging these events. Similarly, 
the v e reaction has threshold of 14.4 MeV and produces two photons from the emitted positron 
annihilation. Also, the 12 B produced in the reaction decays by emitting an electron with a mean 
lifetime of 20.2 ms, providing an additional coincidence signal [23]. With such high thresholds, 
these reactions depend critically on any mixing from the higher energy u x spectra. 

Neutrinos with energies above the 15.11 MeV threshold can excite 12 C into a state that 
decays through emission of a distinctive 15.11 MeV, mono-energetic photon. This signal will 
convey no spectral information about the neutrinos since all interactions produce the same 
mono-energetic photon. This process is neutral current so all flavors contribute equally to the 
cross section. Once again, the high threshold limits the interactions to mainly higher energy 
neutrinos. We have taken the cross section for these reactions from Ref. [28] . 

2. Proton Elastic Scattering: v + p — ^ v + p 

The cross section for this process is fairly large and all flavors contribute almost equally. 
Despite the fact that the scattered protons will undergo a significant amount of quenching due 
to ionization, this process can still produce a large signal for a supernova [27J. Although the 
proton elastic scattering would provide additional statistics, we will only use the 15.11 MeV 
carbon peak for our neutral current signal. This is because extracting the peak from the smooth 
background should be quite easy, while separating the proton elastic scattering from the other 
signals would prove more challenging. 

C. Liquid Argon Detectors 

ICARUS is currently the largest liquid argon detector with a mass of 600 tons, and is 
planned to be scaled up to 3 kton [29]. In ICARUS, ionization electrons from liquid argon are 
detected as they travel along a uniform electric field, providing accurate 3-dimensional event 
tracking and energy reconstruction [30]. One of the most attractive features of liquid argon 
detectors is a strong sensitivity to v e because of the relatively large cross section and low energy 
threshold. An additional benefit is that the charged current and neutral current reactions can 
all be distinguished by gamma emission from the excited nuclei. With electron elastic scattering 
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events producing no coincident photons, all four of the primary reactions can be separated, in 
principle. Two proposed detectors at the 100 kton scale are LANNDD [31] and GLACIER [32J. 

1. Scattering With Argon Nuclei: v e + A0 Ar -> 40 K* + e~ , P e + 40 Ar -> 40 C1* + e + , v+Ax ->Ar* 

Supernova v e and v e both produce charged current reactions in argon with thresholds of 
2.8 MeV and 7.9 MeV, respectively [33]. The distinctive photon peaks produced by the de- 
excitation of the resulting 40 K* or 40 C1* nucleus enables these events to be uniquely identified. 
The first nuclear excited state of argon is quite low, at 1.46 MeV [M], making neutral current 
reactions a particularly strong signal for supernova neutrinos. The subsequent decay of the 
argon nucleus will again provide photon signatures to tag these events, but without a corre- 
sponding electron or positron. The cross sections for charged current events can be found in 
Ref. [33J and neutral current reactions in Ref. [34] . 

2. Electron Elastic Scattering: v + e~ — > v + e~ 

As with previously discussed detector types, electron elastic scattering can also be observed 
in liquid argon detectors such as ICARUS. 

V. EVENT RATES AND RESULTS 

Supernova neutrinos emerge from their neutrinospheres deep in the core of the progenitor star 
and must travel from the extremely dense core to the outer surface of the star. As they travel 
along this density gradient, they undergo oscillations that are dominated by two resonance 
regions: a lower density region (L-resonance) that depends on AmJ 2 and 612, and a higher 
density region (H-resonance) that depends on Am 2 3 and #13. The L-resonance occurs in the 
neutrino sector independently of the mass hierarchy, whereas the H-resonance is significant for 
neutrinos only if the hierarchy is normal and for anti-neutrinos if the hierarchy is inverted. 

The adiabacity of these resonances determines the amount of mixing between the two mass 
eigenstates. In general, this depends on the matter density gradient, the mixing angle, the 
mass difference, and energy of the neutrino. At supernova energies, the L-resonance is perfectly 
adiabatic, leading to essentially full conversion. The H-resonance is also perfectly adiabatic if 
sin 2 2^13 > 10 -3 (referred to as 6*13 large), but will be perfectly non-adiabatic (ie. absent) 
if sin 2 2^13 < 10~ 5 (referred to as #13 small). For intermediate values of #13, the conversion 
probability at the H-resonance is more complicated and depends on energy. For simplicity, we 
will only examine the two boundary cases of perfectly adiabatic and perfectly non-adiabatic for 
the H-resonance. For a detailed discussion of oscillations and flavor conversions in supernova 
neutrinos, see Ref. [35] . 

Once the neutrinos emerge from the remnant star, if there is no intervening exotic physics, 
the mass eigenstates would propagate unchanged to be detected at the earth. We will present 
a method for distinguishing this case from a couple of exotic physics scenarios. 

The number of events from a given reaction can be written as 




(3) 
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TABLE I: The number of events in the normal hierarchy detected in Super Kamiokande for a supernova 
at 10 kpc. The decay scenario is V2, vz—tu\. Here E refers to the kinetic energy of the electron (which 
is the detectable signal) rather than the energy of the neutrino. The Majorana and Dirac rows show 
the number of events for neutrinos of that type if the neutrino masses are hierarchical, producing 
a degraded spectrum for the daughter neutrinos. If the neutrino masses are quasi-degenerate, then 
the daughter spectrum is unaffected (full energy), independent of whether neutrinos are Majorana or 
Dirac particles. The small refers to sin 2 2Qi% < 10~ 5 and large refers to sin 2 2#i3 > 10 -3 . See the text 
for additional initial neutrino parameters. 

with N t the number of targets in the detector, <&(E U ) is the flux of the relevant neutrino at the 
detector, and a(E u ) is the cross section for the process. The limits of the integration are set by 
the relevant range of neutrino energies. We have also assumed 100% detection efficiency and 
neglected the effects of energy resolution. 

A. Results For Water Cerenkov Detectors (Super Kamiokande) 

In Table |TJ we present the numbers of events predicted to be observed from a 10 kpc dis- 
tant supernova in the Super Kamiokande experiment. Here we have assumed a normal mass 
hierarchy and considered four possible exotic physics scenarios. First, we consider the case of 
decays of v% and v-i to v\ without any significant energy losses, denoted by "Decay (Full Ener.)" . 
The columns labeled "Decay (Majorana)" and "Decay (Dirac)" represent the decay scenarios 
in which the neutrino masses are hierarchical, leading to a degraded energy spectrum for the 
daughter neutrinos. If the neutrino masses are quasi-degenerate then the daughter spectrum 
is unaffected (full energy), independent of whether neutrinos are Majorana or Dirac particles. 
Results are shown for both small and large values of #13. In each decay scenario, we assume 
that the propagation time significantly exceeds the (boosted) lifetimes of the decays. Finally, 
the lowest column of the table denotes the case in which all flavor information is lost due to 
quantum decoherence. 

In the upper frame of Fig. [2j we present these results in terms of ratios of event types. The 
first of these, R 2 o, is the ratio of total continuum events (including inverse beta decay, electron 
elastic scattering, and oxygen charged current) above 20 MeV to the total continuum events 
below 20 MeV. The second, -R OX y, is the ratio of total continuum events below 20 MeV to oxygen 
peak events. 

The general features of the plot can be understood as follows. As inverse beta decay events 
dominate the continuum signal, the behavior of R20 is mainly determined by the u e spectrum. 



9 



* 20 




Hierarchical 

Masses 



Normal 
Hierarchy 



Intermediate Masses 



Quasi-degenerate 

Masses 



decays to v 2 

6| 3 small 



10 



No Decay 



1.2 



fl 13 small e 13 large 

1.4 1.6 




2.0 



2.2 



Decoherence 

2.4 2.6 



R 



20 



14 



12 



10 



v 2 decays to y\ 



'y No Decay 

V2 decays to V3 

+ 

013 small 



Inverted 
Hierarchy 



U large 

Quantum ,— j 

+ 
Decoherence 



v 2 decays to V! 



V] decays to v 3 



y? decays to y 3 



I ' 1 Decay- 



i 



y! decays to y 3 



1.0 



1.5 



2.0 



2.5 



3.0 



3.5 



^20 



FIG. 2: Event class ratios at Super Kamiokande for the normal (top) and inverted (bottom) mass 
hierarchies. R oxy is the ratio of total continuum events (including inverse beta decay, electron elastic 
scattering, and oxygen charged current events) below 20 MeV to oxygen peak events. R20 is the ratio 
of total continuum events above 20 MeV to the total continuum events below 20 MeV. The error 
bars represent the la statistical uncertainties. Throughout this paper, black and green error bars 
denote Majorana and Dirac decays, respectively, while red error bars denote full energy decays and 
blue indicate no decay. 



A larger value of R20 indicates that there has been a larger transformation of the initial higher 
energy u x flux to v e through either oscillations in the progenitor star or through exotic physics. 
The oxygen peaks, in contrast, are produced through neutral current reactions of the high 
energy neutrinos. The number of neutral current events will be the same as long the number 
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and spectrum of the daughter neutrinos are the same as for the parent neutrinos. For this 
reason, the largest affect on R oxy comes from scenarios where the decay spectra are degraded. 
In this case, the daughter neutrino only receives a fraction of the parent neutrino's energy, so 
fewer high energy neutrinos exist after the decay. This leads to a large increase in R oxy since 
the number of oxygen peak events decreases significantly. 

At this time, several points pertaining to Fig. [2] should be noted. The quantum decoherence 
signal is well separated from the other scenarios in either hierarchy. In the inverted hierarchy 
we have shown the three decay scenarios that are consistent with the supernova 1987A. For the 
normal hierarchy, only two decay scenarios are presented: ^2,^3 — > v\ and V3 — > V2- Although 
other decay scenarios are possible, those shown predict maximum and minimum values of the 
plotted ratios. Thus any (full energy) decay scenario will be bounded by the four red points 
plotted in the figure, and any degraded decay will fall between those shown on the plot. This 
means that if the neutrinos do decay, this plot should be able to distinguish between the 
quasi-degenerate (full energy decay) and hierarchical (degraded energy decay) neutrino mass 
spectrum. It will be more difficult to discriminate between Majorana and Dirac cases, however. 
Such measurements may even be capable of distinguishing between small and large values of 
#13 at approximately the 2a level. 

If the neutrino mass hierarchy is inverted, the situation is somewhat more complicated. In 
particular, predictions for the ratios shown in Fig. [2] depend significantly on the value of #13. 
Without knowledge of this quantity, most decay and decoherence scenarios we have considered 
cannot be identified using the ratios shown in Fig. [2j The exception to this conclusion are 
decays of the type, v\ — > is 3 , which can potentially be distinguished. 

To overcome this challenge in the case of the inverted hierarchy, we consider signals from 
neutrinos, rather than anti- neutrinos. Unlike anti-neutrinos, the spectrum of neutrinos emerg- 
ing from the remnant star is independent of the value of #13 in the inverted hierarchy. Neutrino 
detection, however, is not the strength of the water Cerenkov detectors. Electron elastic scat- 
tering events, which provide the best measurement of the v e flux, are far less numerous than 
those from inverse beta decay (see Table [T]). They can, however, be separated from the large 
rate of inverse beta decay events by introducing an angular cut. 

In Fig. [3j we show ratios of events chosen to better separate the various scenarios being 
considered in the case of the inverted mass hierarchy. These ratios are essentially the inverse 
of the ratios shown in Fig. |2j but using electron elastic scattering events in place of the total 
continuum events (the reason for inverting the ratio is because the elastic scattering rate is 
quite small, so inverting the ratio gives a larger possible range). The error bars that are shown 
are for a future water Cerenkov detector with 500 kton fiducial mass (rather than the 32 kton 
of Super Kamiokande). 

As can be seen from Fig. [3j these ratios do provide separation between the different scenar- 
ios with the better statistics from a larger detector. Unfortunately, the error bars for Super 
Kamiokande would be approximately four times larger, indicating that Super Kamiokande 
would likely only be able to separate the v\ — > 1/3 decay scenario if the mass hierarchy is 
inverted. It is also interesting to note that for these ratios, a coincidence in the initial spec- 
tra causes the predictions in the no decay and quantum decoherence scenarios to be almost 
identical. The quantum decoherence signal is well separated in both frames of Fig. [2] though, 
indicating that those ratios provide the best opportunity to identify this signal, independent of 
mass hierarchy. 
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FIG. 3: Event class ratios at a future large volume (500 kton) water cerenkov detector in the case 
of the inverted mass hierarchy. i? xy,E is the ratio of oxygen peak events to elastic scattering events 
below 20 MeV. i?20,£ is the ratio of elastic scattering events below 20 MeV to the elastic scattering 
events above 20 MeV. Again, the error bars represent the la statistical uncertainties. 
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TABLE II: The same as Table [I] but for the KamLAND scintillating detector. 



B. Results For Scintillation Detectors (KamLAND) 

Using KamLAND, with a mass of 1 kton, as a representative scintillating detector, we show 
in Table [TT] the numbers of various types of events, in various exotic physics scenarios. Fig. [4j 
we again introduce R20 as the ratio of continuum events (inverse beta decay, electron elastic 
scattering, and charged current events on carbon nuclei) above 20 MeV to continuum events 
below 20 MeV, and compare this to R c , which we define as the ratio of total continuum events 
below 20 MeV to the number of events in the gamma ray peak at 15.11 MeV from the de- 
excitation of the carbon nucleus. 

Again, we find that each of the decay and decoherence scenarios occupies a distinct region 
in the plot. The error bars here a bit larger than in Fig. Rowing to the smaller detector mass of 
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FIG. 4: Event class ratios in the KamLAND detector in the case of the normal mass hierarchy. R c is 
the ratio of total continuum events (includes inverse beta decay, electron elastic scattering, and carbon 
charged current) below 20 MeV to the 15.11 MeV gamma ray peak events from the de-excitation of 
the carbon nucleus. R20 is the ratio of total continuum events above 20 MeV to the total continuum 
events below 20 MeV. Again, the error bars represent the la statistical uncertainties. 

KamLAND compared to Super Kamiokande. The error bars of Fig. [4] would be roughly seven 
times smaller for the planned scintillating detectors in the 50 kton range, however, giving such 
future detectors even better resolving power than Super Kamiokande. 

As in the case of water cerenkov detectors, if the neutrino mass hierarchy is inverted, oscil- 
lations complicate the situation. In fact, even the next generation scintillating detectors in the 
50 kton mass range are not well equipped to distinguish the decay scenarios for the inverted 
mass hierarchy. As discussed in the previous section, neutrino oscillations play the largest role 
in determining the signals in the anti-neutrino sector. Oscillations play a much more minor 
role in the neutrino sector providing the hope that the neutrino flux might provide the ability 
to separate the decay scenarios. 

Unfortunately, the charged- current reactions of v e on carbon nuclei has a very high energy 
threshold of 17.3 MeV, and in the inverted hierarchy, V2 will have the softest spectrum. This 
means that any scenario in which v% decays will have only minor affects on the ratios. We 
also consider the possibility of using electron elastic scattering in the same manner as was 
done for Super Kamiokande. This is not practical, however, because electron elastic scattering 
requires detectors on the order of 500 kton to effectively separate the scenarios (see Fig. |3|, 
approximately 10 times larger than even the next generation scintillating detectors (which are 
of similar size to Super Kamiokande). 

The one scenario that would be very easily identified is again v\ — > 1/3. In this scenario, the 
charged current reactions on carbon would virtually disappear because nearly all of the high 
energy neutrinos will be 1/3, which has a very small projection onto v e . 
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TABLE III: The number of events in the inverted hierarchy detected in the 0.6 kton liquid argon 
detector ICARUS for a supernova at 10 kpc. The decay scenario is vi—±v\. The labels are similar to 
Tables |U and IE 



C. Results For Liquid Argon Detectors (ICARUS) 



Both water Cerenkov and scintillating detectors are much more sensitive to z/ e , providing 
good separation ability between decay scenarios in the normal mass hierarchy. This is because 
the anti-neutrino sector is best for isolating the decay differences from oscillation effects in the 
normal mass hierarchy. The true strength of a liquid argon detector is realized in the inverted 
mass hierarchy. In this case, the neutrino sector provides the best opportunity for isolating 
decay effects. The low energy threshold (2.8 MeV) for charged-current neutrino reactions with 
Argon nuclei gives ICARUS comparatively much better sensitivity to the v e flux over the other 
two detector types. Although the final proposed mass for ICARUS is 3 kton we use the 
current detector mass of 0.6 kton as our model case for liquid argon detectors. 

In principle, ICARUS will be able to isolate each of the four types of interactions in the 
detector. This is another huge benefit for liquid argon over the other types of detectors. 
Some representative results for the number of events in ICARUS for the different scenarios 
in the inverted mass hierarchy are shown in Table III The strong neutrino charged-current and 
neutral current signals are quite noticeable in the table. 

We again introduce two ratios to help to separate the different decay scenarios. We define 
i?30 as the ratio of charged-current neutrino events above 30 MeV to those below 30 MeV. 
We also define R^r as the ratio of total neutral current events to the charged-current neutrino 
events below 30 MeV. These ratios are plotted in Fig. [5] For clarity, the la errors bars for 
the 3 kton mass detector are plotted in the the figure. In comparing Fig. [5] to the lower frame 
of Fig. [2j we find that, in its final form, ICARUS will already be competitive with the next 
generation water Cerenkov detectors in the case of the inverted mass hierarchy. A future liquid 
argon detector of 100 kton would reduce the error bars by about a factor of 6, providing the 
best separation power in the inverted hierarchy by a wide margin. 

The threshold for anti-neutrino charged current reactions on argon nuclei is not particularly 
large at 7.5 MeV. However, the much lower threshold of 1.8 MeV for inverse beta decay ensures 
that both water Cerenkov and scintillating detectors of similar size will have better separation 
sensitivity in the normal mass hierarchy. We have calculated the results in the normal mass 
hierarchy (not shown) for the next generation liquid argon detectors and found similar scenario 
separation to that of Super Kamiokande (see Fig. [2]) . 

All of the results that we have presented thus far are dependent upon the strong hierarchy 
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FIG. 5: Event class ratios in the ICARUS detector in the case of the inverted mass hierarchy. i?3o 
is the ratio of charged-current neutrino events above 30 MeV to below 30 MeV, and Rat is the ratio 
of total neutral current events to the charged-current neutrino events below 30 MeV. The error bars 
represent the la statistical uncertainties, for the planned final detector mass of 3 kton (current mass 
is 0.6 kton). Notice that the quantum decoherence point is coincidentally indistinguishable from the 
no decay scenario. 



of initial average energies within the Livermoore model. We also calculated similar results to 
those presented above using the the softer v x spectrum of the Garching model. Although the 
statistical uncertainties are increased, we still found that the detectors were able to separate 
the scenarios. However, the values of the corresponding ratios between the two models do 
not match well. If a particular model emerges as the accepted one among the community, 
then our calculations could very easily be modified to the accepted model. One way to avoid 
model dependent conclusions is by focusing on the neutronization burst from the supernova as 
discussed in the following section. 



VI. NEUTRINOS FROM THE NEUTRONIZATION BURST 



At the onset of instability, the core begins to collapse, significantly increasing in density. 
Once the density of the inner core reaches that of nuclear matter of order (10 14 g cm" 3 ), the 
core bounces, producing an outgoing shock wave. The shock wave travels towards the outer 
layers of the stellar core, losing energy through the photodissociation of nuclei. The newly 
created free protons have a large electron capture rate, which leads them to neutronize through 
the process e~ +p — > n + u e , resulting in a large flux of v e . These neutrinos are initially trapped 
behind the shock wave, as the density of the core is too high. Once the shock wave has traveled 
outwards to the region of the neutrinosphere (p ~ 10 11 g cm -3 ), however, the neutrinos rapidly 
escape in what is often termed the neutronization burst. 

The presence and structure of the neutronization peak is one of the most robust features 
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of supernova models. In Ref. |36j, the authors show that several different models using many 
different progenitor stars produce very similar peaks. Additionally, the authors examined the 
dependence of the signal in a megaton water Cerenkov detector on several important astro- 
physical inputs. They found that differences in the progenitor mass, three possible models for 
nuclear equation of state for the core, and possible improvements to interaction cross sections 
all led to changes in event numbers that were small compared to the statistical fluctuations. 

The robustness of the neutronization peak across many models provides the opportunity for 
examining exotic physics scenarios without relying on the details of the underlying supernova 
model. One limitation of using the neutronizaton burst peak, however, is that the signal is 
not particularly bright. The neutronization peak typically has a full width half maximum 
of 5-7 ms and a peak luminosity of 3.3 — 3.5 x 10 53 erg s _1 [36]. The total energy released 
is ~3 x 10 51 erg, compared to the much larger ~5 x 10 52 erg of energy in v e for the entire 
supernova. As a consequence, very large volume detectors will be required to discriminate 
between the various exotic physics scenarios being considered here. 

In Ref. p3], the authors examine neutrino decay (for three different decay scenarios) in the 
case of the normal mass hierarchy with hierarchical masses, in which the neutronization peak 
can be identified in the inverse beta decay signal in a future water Cerenkov detector. Similarly, 
we also find that future planned detectors will possess the ability to clearly separate between 
our exotic physics scenarios when focusing only on the neutronization burst. In Fig. [6j we 
plot the same ratios as in the previous sections, but for the neutronization signal only. The 
error bars shown are again statistical and correspond to a 500 kton water Cerenkov detector 
in the normal hierarchy and for a 100 kton liquid argon detector in the inverted hierarchy. An 
important feature of this plot is shown in the normal hierarchy with quasi-degenerate masses. 
If the both V2 and nu^ decay, then the water Cerenkov detector would be able to make a strong 
distinction between the the Majorana and Dirac nature of the neutrinos. 

VII. CONCLUSIONS 

In this article, we have calculated the expected neutrino spectrum and corresponding signals 
from a galactic core-collapse supernova using the Livermoore model. We utilize the extremely 
long baseline to investigate possible signals of various exotic physics models. By introducing 
ratios of event types to distinguish between these various scenarios, we eliminate some pos- 
sible supernova uncertainties (such as total energy released, distance to supernova, etc.). In 
the case of a normal mass hierarchy, Super Kamiokande will be able to distinguish between 
various decay scenarios. In the case of an inverted mass hierarchy, ICARUS is best able to 
discriminate between the possible models. In either mass hierarchy, Super Kamiokande is best 
equipped to identify the signatures of quantum decoherence. We also find that future water 
Cerenkov (~500 kton) and liquid argon (~100 kton) detectors will be able to make strong 
conclusions about these phenomenon in a model independent way by examining the neutrino 
signal associated with a supernova's neutronization burst. 

This work has been supported by the US Department of Energy, including grant DE-FG02- 
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FIG. 6: Event class ratios for the neutronization burst signal. The ratios are the same as in Fig. [2] 
for a 500 kton water Cerenkov detector in the normal hierarchy (top) and as in Fig. [5] for a 100 kton 
liquid argon detector in the inverted hierarchy (bottom). The error bars represent the la statistical 
uncertainties. 

95ER40896, and by NASA grant NAG5-10842. 



[1] J. F. Beacom, N. F. Bell, D. Hooper et at, Phys. Rev. Lett. 90, 181301 (2003). hep-ph/0211305| 



J. F. Beacom, N. F. Bell, D. Hooper et al, Phys. Rev. D68, 093005 (2003). |hep -ph/0307025 



17 



[2] D. Hooper, D. Morgan, E. Winstanley, Phys. Lett. B609, 206-211 (2005). hep -ph/041 0094]; 
D. Hooper, D. Morgan, E. Winstanley, Phys. Rev. D72, 065009 (2005). [hep^ ph/0506091]; 
L. A. Anchordoqui, H. Goldberg, M. C. Gonzalez- Garcia et al, Phys. Rev. D72, 065019 (2005). 
[hep-ph/0506168| ; 

[3] F. Halzen, D. Hooper, Rept. Prog. Phys. 65, 1025-1078 (2002). | astro-ph/0204527l ; T. K. Gaisser, 
F. Halzen and T. Stanev, Phys. Rept. 258, 173 (1995) [Erratum-ibid. 271, 355 (1996)] |arXiv:hep- 
|ph/9410384 |; J. F. Beacom, |arXiv:1004.33Tl"1 [astro-ph.HE]]; A. Burrows, Ann. Rev. Nucl. Part. 
Sci. 40, 181-212 (1990). 

[4] K. Choi, A. Santamaria, Phys. Lett. B267, 504-508 (1991); A. Acker, A. Joshipura, S. Pakvasa, 
Phys. Lett. B285, 371-375 (1992); A. Acker, S. Pakvasa, J. T. Pantaleone, Phys. Rev. D45, 1-4 
(1992). 

[5] K. Nakamura et al. [Particle Data Group], J. Phys. G 37, 075021 (2010). 

[6] J. F. Beacom and N. F. Bell, Phys. Rev. D 65, 113009 (2002) |arXiv:hep-ph/0204111|. 



[7] K. Eguchi et al. [KamLAND Collaboration], Phys. Rev. Lett. 92, 071301 (2004) |arXiv:hep- 
|ex/031004 7|. 

[8] S. Ando, Phys. Rev. D 70, 033004 (2004) l jarXiv:hep-ph/0405200| . 

[9] S. W. Hawking, Commun. Math. Phys. 87, 395 (1982); J. R. Ellis, J. S. Hagelin, D. V. Nanopoulos 
et al, Nucl. Phys. B241, 381 (1984). 
[10] E. Lisi, A. Marrone, D. Montanino, Phys. Rev. Lett. 85, 1166-1169 (2000). 
G. L. Fogli, E. Lisi, A. Marrone et al., Phys. Rev. D67, 093006 (2003). 



hep-ph/0002053 



hep-ph/0303064 



[11 
[12 

[13 
[14 

[is; 

[16 

[17 
[18 
[19 
[20 
[21 
[22 

[23 
[24 
[25 



G. Barenboim, N. E. Mavromatos, JHEP 0501, 034 (2005). |hep-ph/0404014j; G. Baren- 
boim, N. E. Mavromatos, S. Sarkar et al, Nucl. Phys. B758, 90-111 (2006). [hep-ph/0603028] ; 
G. Z. Adunas, E. Rodriguez-Milla, D. V. Ahluwalia et al., Phys. Lett. B485, 215-223 (2000). 
|gr-qc/000602l] ; H. V. Klapdor-Kleingrothaus, H. Pas, U. Sarkar, Eur. Phys. J. A8, 577-580 
(2000). [hep-ph/0004l"23| ; G. L. Fogli, E. Lisi, A. Marrone et al, Phys. Rev. D76, 033006 (2007). 
[arXiv:0704.2568 [hep-ph]]. 



T. Totani, K. Sato, H. E. Dalhed and J. R. Wilson, Astrophys. J. 496, 216 (1998) |arXiv:astro- 
|ph/9710203| . 

S. Horiuchi, J. F. Beacom and E. Dwek, Phys. Rev. D 79, 083013 (2009) |arXiv: 0812.31571 [astro- 
ph]]. 



G. G. Raffelt, M. T. Keil, R. Buras, H. T. Janka and M. Rampp, larXiv:astro -ph/0303226 
R. J. Wilkes, |arXiv:hep-ex/0507097| 

A. Tonazzo [MEMPHYS Collaboration], Nucl. Phys. Proc. Suppl. 168, 363 (2007). 

e.g. , T. Nakaya, presented at the XXth Int. Conf. on Neutrino Physics and Astrophysics, Munich, 

May 2002. 

A. Strumia and F. Vissani, Phys. Lett. B 564, 42 (2003) [arXiv:astro-ph/0302055] . 

J. F. Beacom and M. R. Vagins, Phys. Rev. Lett. 93, 171101 (2004) [arXiv:hep-ph/0309300| . 

A. Kibayashi and f. K. Collaboration, |arXiv:0909.5528| [astro-ph.IM]. 

C. Giunti and C. W. Kim, Fundamentals of Neutrino Physics and Astrophysics, Oxford 2007. 

K. Langanke, P. Vogel and E. Kolbe, Phys. Rev. Lett. 76, 2629 (1996) |arXiv:nucl-th/9511032| . 

E. Kolbe, K. Langanke and P. Vogel, signatures for neutrino Phys. Rev. D 66, 013007 (2002); 

J. F. Beacom and P. Vogel, Phys. Rev. D 58, 053010 (1998) [arXiv:hep-ph/9802424] . 

A. Bandyopadhyay, S. Choubey, S. Goswami and K. Kar, |arXiv:hep- ph/0312315. 

N. Y. Agafonova et al, Astropart. Phys. 27, 254 (2007) |arXiv:hep-p h/0609305 1 . 

T. Marrodan Undagoitia, F. von Feilitzsch, M. Goder-Neff, K. A. Hochmuth, L. Oberauer, 



18 



W. Potzel and M. Wurm, Prog. Part. Nucl. Phys. 57, 283 (2006) [J. Phys. Conf. Ser. 39, 278 
(2006)] |arXiv:hep-ph/0605229| . 

J. G. Learned, S. T. Dye and S. Pakvasa, arXiv:0810.4975 [hep-ex]. 

J. F. Beacom, W. M. Farr and P. Vogel, Phys. Rev. D 66, 033001 (2002) [arXiv:hep-ph/0205220] . 

M. Fukugita, Y. Kohyama and K. Kubodera, Phys. Lett. B 212, 139 (1988). 

A. Badertscher [ICARUS Collaboration], Nucl. Instrum. Meth. A 535, 129 (2004). 

A. Ankowski et al. [ICARUS Collaboration], |arXiv:0812.2373| [hep-ex]. 

D. B. Cline, F. Sergiampietri, J. G. Learned and K. McDonald, Nucl. Instrum. Meth. A 503, 136 
(2003) [arXiv:astro-ph/0105442] . 

A. Rubbia, J. Phys. Conf. Ser. 171, 012020 (2009) [arXiv:0908.1286l [hep-ph]]. 

M. Sajjad Athar and S. K. Singh, Phys. Lett. B 591, 69 (2004). 

I. Gil Botella and A. Rubbia, JCAP 0310, 009 (2003) [arXiv:hep-ph/0307244] . 

A. S. Dighe and A. Y. Smirnov, Phys. Rev. D 62, 033007 (2000) |arXiv:hep-ph/9907423| . 

M. Kachelriess, R. Tomas, R. Buras, H. T. Janka, A. Marek and M. Rampp, Phys. Rev. D 71, 

063003 (2005) [arXiv:astro-ph/0412082| . 



19 



